SU(3) phase states and finite Fourier transform 



Brandon Zanette and Hubert de Guise 

Department of Physics, Lakehead University, Thunder Bay, Canada 

Abstract. We describe the construction of SU(3) phase operators using Fourier-like 
transform on a hexagonal lattice. The advantages and disadvantages of this approach 
are contrasted with other results, in particular with the more traditional approach 
based on polar decomposition of operators. 



1. Introduction: complementarity and the Fourier transform 

The idea of complementarity in quantum mechanics goes back to Bohr and his attempt 
to explain wave-particle duality. The concept was sharpened by Pascual Jordan, who 
stated pQ that 

For a given value of x all values of p are equally possible. 

This formulation automatically singles out the Fourier transform connecting operators 
like x and p as their respective (generalized) eigenstates satisfy 

(x\p) ~ e ixp/h \(x\p)\ 2 = constant (1) 

The concept is not limited to continuous systems but also exists in finite dimensions. In 
this contribution, we will discuss the construction of SU (3) phase operators, which are 
expected to be complementary to number operators. This contribution emphasizes the 
importance of the finite Fourier transform, and in particular a new type of generalization 
of the Fourier transform that is constructed to preserve the symmetry of a hexagonal 
lattice, which is the natural (discrete) lattice to describe states appropriate to the 
description of a collection of three-level systems. Our approach should be contrasted 
with the approach of Dirac [2], which emphasizes polar decompositions, and which has 
been applied to SU{2) and other systems in [5] and [3]. 

2. Two examples 

Consider first a spin-| system, taking as operators the Pauli matrices a x , a y and a z . The 
eigenstates {|+)z, \—)z} OI °~z and the eigenstates |— ) x } of o~ x are complementary: 

\x{+ I +)z\ 2 = \ x (- I +}z\ 2 = \x(+\ ~)z\ 2 = \ x (- I -)z? = \ = constant .(2) 
The eigenstates of a z and a x are related by a finite Fourier transform: 
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The operators o~ x and o z are said to be complementary. The same property holds for 
the pair o~ y and a z and for the pair o y and a x . The transformation matrix connecting 
any two sets of eigenstates remains a finite Fourier transform. 

A similar construction exists for a three-level system (or qutrit). Defining 
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w = e lm,i , (4) 

and writing their respective eigenstates as {\0 Z ), |1 2 ), |2 2 )} and {|0 X ), \2 X )} we find 
for instance |(la;|0 2 )| 2 = |(2 a .|2 z )| 2 = | with all other such overlaps constant. Here 
again, the eigenstates of X and Z are related by a finite Fourier transform: 
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This is a good point to mention some of the properties of the finite Fourier matrix 
F. It is unitary, which implies 

»U (6) 

i 

(This would be orthogonality under integration in the continuous case.) F 4 = t, and 
its entries are characters of finite Abelian groups. In dimension n: 

Finally, have constant magnitude, connecting with Jordan's definition of 

complementarity. 

3. SU{2) phase states 

Following Dirac (2] and others [3], phase operators in su{2) (and other) systems are 
constructed by writing the matrix for S + (or S^) in polar form, viz. 

3-1 

S+^^c m \j,m+l){jM=E-D (8) 

where D is diagonal and E is a "phase" part, containing entries which produce the 
shifting action of S + on the basis states. The operator E is expected on physical 
grounds to be complementary to the diagonal operator S z . 

Geometrically, the set of eigenvalues {m} of S z acting on number states \jm) are 
equidistant points on a line and the action of the ladder operators S± takes a point m 
to its neighbor mil. The action of S+ is pictorially represented in Fig.. 

Because \jj) is killed by S + , the rank of S + is one less than the dimension of the 
system, so that E is not completely defined: we can adjust the entries in one line. The 
usual choice makes E cyclic (i? 2j+1 = 1) so it generates an Abelian group of order 2j + 1: 
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Figure 1. The shift action of the su(2) raising operator 



This E is unitary, and can be written in the form E = e 1 ^, with (p the putative hermitian 
phase operator. The eigenvectors of E are eigenvectors of ip and defined to be the SU(2) 
phase states. The components of the m'th eigenvector \ip m ) are just elements of a Fourier 
matrix F. Thus, the phase eigenstate \<p m ) is given by 



i) = ^2F mk \jk} 



£irikm/(2j+l) 



(10) 



4. 377(3) and SU(3) phase states 
4-1. Geometry of SU (3) states 

The algebra su(3) appears naturally in the construction of number-preserving transition 
operators for three-level systems. There are six transition operators, usually denoted 
by Cij = a\aj for i ^ j = 1,2,3, and two population differences hi and 
h,2 = a\a2 — a\a%. The states 1 200) and 1 110) , for instance, correspond to the pairs (2, 0) 
and (0, 1) of population differences. Pairs are located on a hexagonal lattice with basis 
vectors u\ and U2 as illustrated on the left of Fig|2j 





Figure 2. Left: The population differences (rii — n 2 , n 2 — n 3 ) for the states |7T. 1 n 2 n 3 ) 
with rii + n 2 + n 3 = 2. (a, b) is located at auj\ + buj 2 with w 1; w 2 the lattices vectors. 
Right: The graphical representation of the shift action of the ladder generators of su(3) 
on the hexagonal lattice. The ringed dot at the center represents the two diagonal 
population difference operators hi and h 2 , which do not shift the basis states. 



The action of on lattice points is illustrated on the right of Fig. |2j Basis vectors 
cti and a 2 associated to the operators C*i 2 and C 2 3 respectively are dual (reciprocal) to 
the lattice vectors uj\ and U2, respectively, as illustrated. Using the hexagonal geometry, 
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two points (a, b) and (c, d) corresponding to two pairs of population differences differ by 
an integer combination of the vectors a\ and a 2 - The action of on the state \n1Ti2Ti3) 
is to translate the point (jii — n,2,n2 — TI3) by the vector associated to to the point 



(n\ — n' 2 , n' 2 — n'i), so that, for instance 



"3 J J 

C 12 |110> 



|200) (0,1) 1 y (2,0) =a x + (0,1). 



The central ringed dot represents the two diagonal population difference operators hi 
and h 2 - There should be one phase operator conjugate to each hi. 



4-2. Two solutions; boundaries 

If we approach the construction of SU (3) phase states using polar decompositions, we 
are faced with an interesting problem. Because there are two basic shift directions, ot\ 
and CK2) each one of C12 and C23 will come with its own set of not necessarily mutually 
compatible boundary conditions. 

In the simplest case of the states 1 100) , |010) and |001), the shift matrices £12 and 
E23 that enter in the decompositions of C12 and C23 respectively contain two lines that 
cannot be uniquely determined. These matrices can be completed in two different ways. 
First, we can write 
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[|010)(100| + |100)(010|] + |001)(001| 



|100)(100| + [|010)(001| + |001> <010|] 



(12) 
(13) 





Figure 3. A graphical representation of how cyclic boundary conditions can be 



imposed to complete the matrix E\ 2 to yield Eqn.(12| 



This kind of solution also exists for more general cases where n\ + n 2 + 713 > 1. 
It consists in considering subsets of states with the same value for n 3 - such subsets of 
states fall on lines parallel to the o>i direction - and following the procedure of SU (2) on 
each lineto obtain E12. Similarly, by considering subsets of states with the same value 
of rii - these states now lie on lines parallel to the «2 direction - we can follow the 577(2) 
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procedure for each line and obtain E23. However, one feature of this solution, already 
present in Eqns.(12) and (13), is that the unitary phase operators do not commute: 

[E12, E 2 s\ 7^ 



j,0<p 



i(e^i2+7<^23) 



(14) 



This in turns implies that the phases are not additive. 

For the case of the states 1 100) , |010) and 1 00 1) , it is possible to find shift matrices 
E\2 and E23 compatible with the polar decomposition of the respective operators and 
so that [E12, E 23 ] = 0. These matrices are 
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Figure 4. A graphical representation of how cyclic boundary conditions can be 



imposed to complete the matrices E\i and E23 to yield Eqn.(15) 



Note that Ei 2 is just the operator X of Eqn.(|4j) while E 23 = X 2 . Clearly, E\i and E 23 
commute. However, we have not been able to find similar solutions for sets of states 
with rii + n 2 + n 3 > 1. 

4-3. Finite Fourier transform on a hexagonal lattice 

As an alternative to the construction based on polar decomposition, we look for a finite 
Fourier transform adapted to the discrete hexagonal symmetry natural to SU (3) states. 
Such an FFT was proposed in [5] and will be adapted to our needs. 

We start with the physical states \nin 2 n 3 ). The procedure of [5] requires that the 
"data points" be in the first hextant of the lattice, so we find a rigid displacement of the 
set of population differences {n\ — n2,n 2 — n 3 ) corresponding to the physical states so 
that every pair {n\ — n 2 , n 2 — n 3 ) is mapped to a single point in the first hextant. One 
can show that the rigid displacement is a linear transformation comprising a translation, 
a rotation and a change of scale of the original pairs of points. The final result of the 
sequence is 

\nin 2 n 3 ) h)> (m - n 2 , n 2 - n 3 ) h)> (m, n 2 ) (16) 

An example of result is given in Fig|5] We obtain for each point (a, b) in the first extant 
its orbit, i.e. the set of points obtained by considering reflections of (a, b) through 
mirrors perpendicular to a± and a 2 . Depending on the value of a and b, an orbit may 
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Figure 5. An example of how the pairs of populations differences (p, q) obtained from 
the states \n1n2n3} with n\ + n 2 + n 3 are mapped to the first hextant. 



contain 1, 3 or 6 points. The orbit for the points (2, 0) and (1, 1) are illustrated in Figj6j 
Each orbit is labelled by its starting point (a, b) in the first extant. There is the same 



\ /H 


.2/ 














)/-2,1)V 










(\ 


■i) 


V' 7 \ 



Figure 6. The orbits of the points (2,0) (in red) and (1, 1) (in green). 

number of orbits as the number of states. Each orbit is used to construct a so-called 
orbit function 

X(ab)(ni,n 2 ) ~ LO (2a+b)n 1 +{a+2b)n 2 UJ {b-a)n 1 + (a+2b)n 2 _|_ w (2a+6)ni+(a-6)n 2 

_j_ ^-(a-b)ni-(b+2a)n 2 _|_ ^- (a+2&)m+(6— a)ri2 _|_ ^-(2b+a)ni-(f>+2a)n 2 (17) 

with (n 1; n 2 ) points in the first hextant. The functions \ as closely related to characters 
of elements of finite order of SU(3). 

It is essential to rigidly translate the population differences of physical states. Two 
states 1^1^2^3)5 \ n i n 2 n 3) which differ only by a permutation of n\, n 2 , % yield population 
differences (rii — n 2 , n 2 — n%) and (n^ — n 2 , n 2 — n' 3 ) that are on the same orbit and so 
produce identical functions x- It is only once the population differences have been 
translated to the first hextant that (ni,n 2 ) and (n'^r^) will lie on different orbits. 

The functions x need to be properly normalized and weighted as described in [5] 
but, once this is done, they satisfy an orthogonality relation 

^ (X(a,b)(ni,n 2 ))* X(a',b>)( n l' n 2) ~ $aa>0~bb> ■ (18) 

ni,n 2 

The orbit functions can then be used to obtain a Fourier matrix 

/ Xa 1 ,6 1 (si,S 2 )l ... Xa 1 ,6 1 (si,S 2 )fc' 

- ; ••. ; I no; 

\Xa k ,b k (s U S2)l ••• Xa k ,b k (Sl,S 2 ) k , 
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So defined the matrix F immediately satisfies the majority of the conditions given at 
the end of Sec. 2. In particular, for the set of states { 1 100) , 1 010) , |001), the matrix F is 
exactly that of Eqn.(j5j). However, for other states with n\ + n 2 + n^ > 1, the matrix F no 
longer contains entries of the same magnitude. For instance, using the states \nin 2 n^) 
with Tii + n 2 + n 3 = 2, we find 
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(20) 



57/(3) phase states 

Now define SU(3) phase states as transforms of the shifted population difference 
eigenstates: 

I ?7l j *72)(n 1 ,n 2 ) = y^^(m,n 2 ),(ti,ta) 1*1) fr?) (21) 



*1,*2 



Phase "operators" are conjugate to population difference operators: 

771 = Fhi F~ l , 772 = Fh 2 F- 1 . 
Since [/ii, hi] = 0, we recover [771,772] — 0: phases commute. 



(22) 



4-5. Complementarity and number difference distributions 

To get insight into what phase states "look like" we consider the probabilities 
\(ni,n 2 I rji, r] 2 )\ 2 . Recall the correspondences \nin 2 n 3 ) -H- (ni,n 2 ) between physical 
states and their translated population difference. Thus 

|100)o(l,0) |010) o(0,l) |001) o(0,0). (23) 

For these states we find, for every {n\,n 2 ) and every (a,b): 

\((n h n 2 )\[F\(a,b))]\ 2 = l (24) 

This is no surprise as \F(mn2),(a,b)\ 2 = 3 f° r this case. 

For ni+n 2 +n 3 = 2, it is convenient to construct probability histograms for points in 
the first hextant. Using as input state any one of the states 1 200) -H- (2, 0), 1 020) -H- (0, 2) 
or 1 002) <r» (0, 0), we find two and only two possible amplitudes, as illustrated with two 
different colors on the left of Fig|7| The corresponding histogram for any one of the 
input states |110) O (1,1) or 1 101) o (1,0) or |011) O (0,1) is on the right of Fig(7| 

For a given input state not every Fourier component has the same amplitude: 
complementarity in the sense of Jordan is lost - as expected since |-F(mn 2 ),(a,fc)| 2 is no 
longer constant. However, points (ni,n 2 ) and (n'i,n 2 ) with equal amplitudes in the 
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Figure 7. Left: probability histogram for any one of the input states 1 200) , 1 020) or 
1 002) . Right: probability histogram for any one of the input states 1 110) , 1 101) or 1 011) . 
Columns of the same color have the same height. 

first hextant correspond to physical states |n 1 n 2 n 3 ) and In^n^n^) which differ by a 
permutation of their entries. 

For a generic input state, like |0, 21, 9), the probability landscape is rugged without 
any special features. However, for input states of the type |iV00) or |0iV0) or |00iV), 
which are mapped to the corner edges of the first hextant, we find that the probability 
landscape is remarkably regular. The probability landscape is asymptotically flat, 




Figure 8. Probability landscapes for input states |5, 0, 0), 1 10, 0, 0) and 1 20, 0, 0). 



meaning that, in the large N limit, the phase states F\N00) etc are asymptotically 
conjugate to the Fock states \nin 2 n 3 ). 

4-6. su{3) phase operators 



The phase operators 771,772 of Eqn.(22) generally have "complicated" expressions. In 
spite of this, we have found the following observation to hold. If we evaluate the variances 
A771 and A?72 using the physical states \n1n2n3), the smallest variances always occur for 
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the states |iVOO), 1 A r 0) or |00iV). The landscape of variances of fji evaluated in the 
physical states |n 1 ra 2 n 3 ) is illustrated in Figj9]for n\ + n 2 + n 3 = 5 and 30. 




Figure 9. Variance of the phase operator r)i calculated using the physical states 
\n\n2nz). Left: ri\ + ri2 + 713 = 5. Right: n\ + ri2 + n% = 30. 



5. Conclusions 

The polar decomposition of operators in SU(3) produces phase operators that are 
ambiguous and not unique: in general, non-commuting raising operators lead to a 
decomposition that produces non-commuting phase operators. Moreover, this approach 
produces an "exponential phase" rather than a phase operator. 

We can obtain hermitian commuting "phase-like" operators by using of symmetry- 
adapted FFT. The procedure is mathematically systematic but not very intuitive, 
and we loose the connection with complementarity. With this approach the physical 
states |iV00), |0iV0) and |00iV) stand out as having unexpected properties of asymptotic 
complementarity. The variances of the phase operators evaluated using those states is 
always the smallest. 

An unanswered question (not discussed in this contribution) is the difficulty in 
imposing correct cyclic boundary conditions on the phase operators themselves once 
they are exponentiated. 
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